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In this paper we study lepton flavor violating semileptonic decays of heavy baryons in a framework 
of the model mm with large extra dimensions and a single generation . Resulting branching ratios 
for these decays mediated by Kaluza-Klein modes of photon and Z-boson are presented in a model- 
independent form. 


I. INTRODUCTION 
A. Preliminaries 

Standard Model (SM) provides extremely accurate de¬ 
scription of the particle physics. Therefore, it can be 
regarded, at least, as an effective low energy theory. To 
probe physics at higher energies, one can study indirect 
effects caused by yet undiscovered heavy particles pre¬ 
dicted by various models. 

Hereafter we will concentrate on the particular model, 
suggested in the series of works mm, which employs 
the concept of large extra dimensions [3] to solve some 
problems arising within the SM, such as generation repli¬ 
cation, fermion mass hierarchy, and gauge hierarchy. The 
model introduces two additional spatial dimensions form¬ 
ing a compact manifold (sphere). The radius of the 
sphere, R , can be understood as an effective scale of 
gauge fields localization [2]. 6-dimensional fermions, on 
the other hand, are localized in a four-dimensional core of 
an Abrikosov-Nielsen-Olesen vortex [5], which is formed 
by two auxiliary fields. Three zero modes of the single 
fermion generation in six dimensions give rise to three 
four-dimensional generations of the SM fermions [Ij. 

This model can be directly probed, since it implies the 
existence of the 4-dimensional fermion currents that are 
forbidden in the SM. This fact was used in Refs. mi, 
where the possible influence on K-, D- and B-meson de¬ 
cay physics was investigated. The experimental bounds 
on the branching ratios of such processes lead to the con¬ 
straints on the parameters of the theory, namely, on the 
typical size R of the localized gauge zero modes. The 
strongest bound found so far, which is 

7 > 64 TeV, (1) 

arises from the non-observation of K° —> pe decay |4j. 

Since the vast field of baryon physics remained un¬ 
touched in the previous works, here we concentrate on 
the processes involving heavy baryons. The goal of the 
present paper is to calculate the partial decay widths of 
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forbidden baryon decays. In the framework of model with 
large extra dimension discussed above, the constrains on 
the decay widths can be translated to the bounds on Rr l . 

Analytical expressions for semileptonic baryon decay 
rates were calculated recently, see, e.g. Ref. [7]. In this 
work we take into account both weak and electromagnetic 
channels which makes the calculation more complicated. 
To obtain squared amplitudes, we use FeynCalc package 
0 . 

B. Model features 

Here we describe a number of model features that are 
crucial for the following calculations. Detailed descrip¬ 
tion and derivation of these features can be found in Refs. 

mm®. 

The main feature is the existence of heavy Kaluza- 
Klein modes of the SM gauge fields. Indeed, due to pe¬ 
riodical boundary conditions on the sphere of extra di¬ 
mensions, 4-dimensional gauge field forms Kalutza-Klein 
tower of heavy excitations. Zero mode corresponds to the 
SM gauge held, and other modes are massive and sepa¬ 
rated from the lowest mode by the mass gap M ~ 1 /R. 

Higher excitations of the SM gauge fieldscan provide 
horizontal transitions between different generations of 
fermions [1] . These transitions can be observed in Lepton 
Flavor Violating (LFV) processes as well as in processes 
with Flavor Changing Neutral Currents (FCNC). 

Interaction of the fermion current j u ' with the higher 
excitations of the SM gauge boson is given by 4- 
dimensional effective Lagrangian ® 

C eff = g ■ Tr(A%), (2) 

where g is an appropriate gauge coupling constant (i.e., 
e for photon). Note that A M and ,j /( are matrices with 
indices m, n enumerating generation number: 

au pi A au Jr ; 1 ’ 2 au \ 
- C, 11-'N,0 ^ 13 ^ 1,2 \ 

E 1 £a^ Ej$A? t0 E 1 £a \, ( 3 ) 

14X2 E l i A% E l £A» 0 ) 

and 

q A 4 — /it n 
Jmn ~ Uj m u 
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where a n are two-component Weyl spinors and E m 
are overlap factors calculated in BNBi- Index l in @ 
enumerates modes A/ with four-dimensional masses 
m 2 = /(/+ 1)^-2. 

The situation becomes more complicated, if one also 
takes into account quark mixing. For our consideration 
it is important that the amplitudes of the reactions that 
change the total generation number AN, are addition¬ 
ally suppressed compared with those which keep it un¬ 
changed, see Ref.0 for details. 


Effective four-dimensional Lagrangian 

In what follows, we use the effective four-dimensional 
FCNC Lagrangian in the form 

^eff = (jZ^mqn + 

GwZ^ n ( (4) 

where A l f nn and Zl^ n stand for massive neutral gauge 
bosons, Q em and Gw are corresponding coupling con¬ 
stants; j e ^ mqrh ,j*J mln and j^q mqn ,j^ Tmln are electro¬ 
magnetic and weak currents which are supposed to have 
the same form as the allowed within the SM currents. 
Namely, for the neutral weak currents we have, 

J% mln = (“ + i) ImOjiln + l;lmO?ln + h.C., 
jZi m q r , = ( T 3 - €Q) VmO^qn - £,Qq m 0^q n + h.c., 


O 


R,L 


1 ± 75 
= 2 


where £ = sin 2 9w, Tz marks the component of the quark 
doublet and Q is the charge of q m in |e| unit. Electro¬ 
magnetic currents read similarly, 


Cu = ! ">^ + kl 

3^ m q„ = Qlm'lu-Qn + h.C. 


( 6 ) 


Note that there is no summation over generation indices 
m,n in Q-Q. 

Lagrangian Q has a rather general form which makes 
our results applicable to any model with FCNC and LFV 
processes by the appropriate choice of the coupling con¬ 
stants G e and Gw- In the model under investigation, 
the matrix determines how heavy gauge bosons are 
coupled to the FCNC. From Eqs.(pl) and pj, we have in 
Eq.@, 


GemAmn — e ^ 


GwZmn = 


1=0 


9 


(7) 


2 cos Ow 


E rpl,l yl.l 


1=0 


II. BARYON DECAYS 
A. General consideration 


We are interested in semileptonic baryon decays. Ac¬ 
cording to the bounds obtained in Refs. Hi, the inter¬ 
mediate bosons are heavier than 64 TeV. This implies 
that one can neglect the momentum carried by these 
bosons, and turn to effective four-fermion interaction. 
The amplitude of B t —» Bf£ m £ n transition then reads, 


M = G em Hi/ ■ j *7 , + G W H»_ 




V—A ' iZuG 


( 8 ) 


with G em = Gem/M 2 and Gw = Gw/M 2 being the ef¬ 
fective four-fermion coupling constants, where the mass 
scale M may be the same or different for electromagnetic 
and weak channels. For our model, one should also take 
into account the infinite number of higher excitations, 
that will be done in the next sections. The hadronic 
current in Eq.Q can be written in the following form 

mm), 


K_ a = B fin (F^(q 2 ) + 75 G^ 2 )) B in , (9) 


where 


( q 2 ) 

(q 2 ) 


/ifaV 

9i(q 2 b M 


/ 2 (g 2 ) 

m B in 


<J^q v + 


m Bin 


Hq 2 ) 

-1 

m B in 

gojg 2 ) 

m B irl 




( 10 ) 


Here q is the momentum transferred to leptons and the 
dimensionless quantities /,;, g t are the hadron form fac¬ 
tors. To the leading order in \q\/ms, we have fi(q 2 ) ~ 
/i(0), 9i{q 2 ) ~ 3i(0), and /* « 0, g t « 0, i = 2,3. 
The next order approximation is fi{q 2 ) « fi{ 0), gi(q 2 ) ~ 
gff 0), * = 2,3, and [5! 


fi(q 2 )* 

gi{q 2 ) ~ 


/i(0) 

(1 + q 2 / m y) 2 
gi(0) 

(1 + q 2 /m 2 A ) 2 


( 11 ) 


Constants m A ,my as far as the values of the form fac¬ 
tors at zero momentum must be specify individually. Ex¬ 
pression for Hy can be obtained from Eq.Q by setting 
G„ = 0. 

For the decay width we have, as usual ms, 

dT= {ly 32M|. \WdmUrr4s, (12) 


where M Pjj is the mass of decaying particle, and m 2 j = 
(Pi +pj) 2 , see definitions of final-state momenta in sub¬ 
section |HB] As was already pointed out, the amplitudes 
with AN ^ 0 are suppressed. We will, therefore, concen¬ 
trate processes with AN = 0. Furthermore, we restrict 
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ourselves with the processes carried by the neutral cur¬ 
rents described above, since we wish avoid neutrinos in a 
final state and detect safely all products of the reaction 
in order to be capable to set constraints on the flavor 
violation. 

As we have already explained, all calculations provided 
below do not depend on the details of the theory whose 
low-energy limit we investigate, so, with obvious replace¬ 
ments of the coupling constants, the decay rates of the 
baryons will be the same in any other theory allowing 
such transitions. 


B. An example: decays of S-baryons 


As an example, let us perform the detailed calcula¬ 
tions of the semileptonic decays of S-baryons. The triplet 
of strange E-baryons is (E + , E _ , E°). For E + and E° 
there are, e.g., reactions E + —► pe + and E° —> ne + , 
which both go through the neutral currents. Two other 
similar reactions E + -A pe~ g + and E° -A ne~ p + have 
|AAf| = 2, and, hence, their decay rates are suppressed 
additionally by a factor ~ 10 -4 , so, we do not consider 
them. Reactions with r-leptons in the final state are im¬ 
possible due to kinematics. 


E + -a pe + g decay 


Both the neutral boson modes Z A and photon modes 
A^'-j contribute to the amplitude of E + -a pe + decay. 
The part of the effective Lagrangian Q, corresponding 


to this process, is 

OO 

£'St = ‘ E^ud-ife+iS,) 


1=0 


9 


2 cos 9 W 


OO 

SdAbfe+ib.). (is) 

(=0 


where and j™ ls , j^ e are the appropriate 

electromagnetic and weak currents. From Eqs. © and 
([b]), we have, 




•CL = - 3..*. = ( -7 + o£ + 775 ) s, 


iW 


1 1 


1 


4 3 4 

j%pe = ~~Pl^ j^-pe = ( - 7 + £ + 775 ) e - 


In the four fermion approximation, the amplitude 
reads, 


V • J u 1. 




+ G W Hy_ A ■ j W 




(14) 


Since p and E + belong to the same SU(3)-octet, we can 
apply the group approach to estimate the relevant form 
factors in Eqs.(10l to the leading order in \q\/m-z jT2] , 
This gives /i 1.2. 

In order to account the infinite tower of heavy exci¬ 
tations of the gauge bosons, we determine the coupling 
constants in Eq. (14) in the following way: 


G w = 


G = e 2 V 

u em ° / J 


2 cos 6w 

(dl) 


s 00 f E la 1 ) 2 

Y V V = V2G F m 2 z (R 2 , 

z —^ mr 


1=0 

2 


= e 2 (R 2 , 


1=0 


where ( represents the sum of the series [4] 

K) r 


(15) 


c-E 


^ 10 + 1) 


0.4. 


We define momenta of initial and final particles in the 
following way: E + (P) -a P + {pi)p~ i,P2)e+[po)- Squared 
amplitude, averaged over spins of the initial baryon and 
summed over spins of the final particles (14), is: 


\M\ 2 = 2 [m p Mzp 2 ■ P 3 (4/i (1 ^ 4 OG em G w - 8 f 2 G 2 em - (8£ 2 - 4£ +1) (. f 2 - g 2 ) G 2 W )) + 

(4 hG em G w (A(4£ - 1) + 9l ) + 8 f 2 G 2 em + G 2 w (2f l9l (4£ - 1) + f 2 (8{ 2 - 4£ + l) + g\ (8{ 2 - + l))) + 

Pi ■ P 2 P ■ p 3 (2figiG w (—2 G em — 4 £Gw + Gw) + fi (4(4£ — l)G em &V+ 

8GL + (8£ 2 - 4£ + l) G 2 w ) + g\ (8>e - 4£ + l) G 2 W )] . (16) 
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We have neglected electron mass in the last equation. In 
order to compute the decay width, we rewrite Eq. (16) 
in terms of m\ 2 and TO 23 , and, using Eq. (12), we obtain 
numerically: 


r ~ t w g 2 w + r w , em GwG em + r em G 2 m , (17) 

with 

T w ~ 6.9 • 10“ s GeV 5 , 
r w ,em^-1.5-10- 7 GeV 5 , 
r em ~ 3.2 • 10“ 7 GeV 5 . 

For the value of 1/R consistent with the restriction ([!]), 
one can find for the branching ratio, 

Br{ E+ -»■ p+p~e + )\ r< r k ^ < 2.4 • 10“ 15 . (18) 


£° —> ne + fj, decay 


This process is also based on the reactions s —1 d + A^ 7 
s —} d T so it is calculated very similarly to E + —> 
pe + . The only difference is in the appropriate form 
factors, for which we have, following to [T2] : /i a: 0.7, 
gi ss 0.2. Finally we get: 

T w ~ 1.2 • 10“ s GeV 5 , 
r w ,em — —2.2 • 10“ s GeV 5 , 
r em ~ 1.6 • 10“ 7 GeV 5 . 

For the value of 1/i? consistent with the restriction ([I]), 
one can find for the branching ratio: 

Br(E° np-e + )\ R<RK ^ e < 1.3 • 10“ 24 . (19) 


The reason for the sharp difference between the low limit 
restrictions (181 and (19) is in the distinction between 
mean lifetimes of the baryons: t^+ ss 0.8 • 10 “ 10 s and 
ty ,o ~ 7.4 • 10“ 2 O s. Due to this difference, the reaction 
E + —> pe + g- is more preferable for searches in experi¬ 
ments. 


C. Charmed and beauty baryons 


Charmed baryons 


Among the charmed baryons, the best restrictions on 
the forbidden semileptonic decays so far have been ob¬ 
tained for A+-baryon. To avoid the additional suppres¬ 
sion, we will consider e —> p lepton current transition, 
since c —> u transition between the first and the second 
generations of quarks have AiV = 1. The interaction 
vertex for these reactions has the following form, 


= +j ™) + 


1=0 


9 


2 cos 0w 


J2 E ti Zl l l 1 ti™uc + 3%j I eJ, 


( 20 ) 


1=0 


with 


2 - 


= ^ 7 ^, j W c = ^ ~ - ^75 ) c, 

1 


A era _ 


-P7v e > j^jie = M7/X ( ~ + £ + ^75 


This leads to the following transitions: 

A+^p, E+(2455)- 
2+ E+, 2° -► E° 


P, ^(2455) 

-7 2 °, 


( 21 ) 

( 22 ) 


with the amplitudes given by the expression (14). Due to 


the significant difference between masses of the initial and 
final baryon states, we cannot use the approximation for 
the form factors j)(0) « 0, <p(0) ~ 0, i = 1,2. We then 
use the values obtained in Ref. [10] under the condition 
that the decay of the baryon is carried by the “heavy 
quark” —►“light quark” transition. The coefficients m Al 
my for the charmed baryons in the dipole approximation 
( [TT]) of /i, gi are equal to [9] 

m v = 2.11GeV, rn A = 2.54GeV. 

We perform calculations of the squared amplitudes with 
the help of FeynCalc package [ 8 ] . The resulting expres¬ 
sion consists of 1042 terms, so we do not quote it here. 
Finally, we perform numerical integration over the phase 
space. Results for the decays listed in Eq. (21) are sum¬ 
marized in the Table |U 


TABLE I. Decay widths and branching ratios. Lw, r Gm ,w 
and r em are defined in Eq. ( 17 1 


Decay 

r w ,GeV° 

rw,em,GeV 5 

E cm, 

GeV b 

Br\ R<RK ^ ue 

A+ pge + 

7 . 7 - 

10“ 4 

7 . 2 - 

10““ 

5 . 3 - 

10““ 

4 . 6 - 

10 44 

E+ ->ppe + 

5 . 8 - 

KT 15 

2 . 5 - 

io- H 

8 . 7 - 

10““ 

4 . 4 - 

10“ 2 “ 

E“ —¥ nge + 

1 . 4 - 

10““ 

4 . 4 - 

io-“ 

2.0- 

10““ 

2.0- 

10“ 24 

E+ -> T, + ge + 

3 . 7 - 

10““ 

7 . 7 - 

10~ b 

1.1 ■ 

10““ 

5 . 7 - 

10“ ib 

Sc -»■ E >e + 

1.8- 

10““ 

3 . 8 - 

10'“ 

5 . 7 - 

10““ 

2.8- 

10“ 15 

fie ->■ H >e + 

1.0- 

nr 5- 

1.1 ■ 

10 & 

3 . 0 - 

10““ 

5 . 8 - 

10“ iY 

A 1 ,) S> E >e + 

9 . 3 - 

Ttr 4 ” 

3 . 0 - 

10~“ 

6 . 7 - 

10““ 

2.34 

10“ 12 

A“ —> nre + 

1.0- 

10““ 

3 . 3 - 

10'“ 

7 . 6 - 

10““ 

2.6- 

10“ 14 

v + 

^ fie 

1.0- 

10'“ 

3 . 3 - 

io-“ 

7 . 4 - 

10““ 

2 . 7 - 

10“ 12 

E5 -► E u re+ 

5.1 • 

10~ 4 

1.6- 

10~“ 

3 . 7 - 

10““ 

1 . 3 - 

10“ 12 

Eft” - 5 - E“/re + 

5 . 7 - 

10'“ 

7 . 2 - 

10'“ 

1.6- 

10“ 2 

7 . 3 - 

10“ 12 

—y fie + 

6.8- 

10“ 

2.2 • 

10“ 2 

5.1 • 

10“ 2 

1.2 • 

10“ 14 

DjC —¥ E“re + 

8.1 • 

10“ 

2 . 7 - 

10“ 2 

6.2- 

10“ 2 

1 . 5 - 

10“ 41 


Beauty baryons 

Forbidden decays of beauty baryons occur by b —> d 
and b —> s transitions. Their consideration repeats the 
previous discussion. We are interested in e —1 r lepton 
current transition for b —> d reaction, and e —> g transi¬ 
tion for b —> s one, which are | AiV| = 0. The appropriate 
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parts of the effective Lagrangian are similar to (13), (20) 


with the obvious replacements. The constant ( in the 
expressions for the coupling constants (15) turns to [5] 


for b —> d, e 


™ ( E ls) 

t^ + 1 2 3 4 5 6 ) 

r transitions and 


0.27 


SM but that, however, can take place in its certain exten¬ 
sions. Since the processes changing the total generation 
number AN are suppressed in the model, we considered 
only those with AN = 0. The restrictions on 1/R, cal¬ 
culated in this work, appeared to be weaker than the 
similar restrictions obtained from B- and, especially, K- 
meson decays [3]. The current experiments, including 
LHCb, could not provide sufficient statistics to increase 
the low limit on i? -1 up to the level ~ lOTeV. 


c" = E 




E 1 ’ 1 E 1 

■ c '12 ^23 

1(1 + 1 ) 


0.47 


for b —> s, e —>■ n transitions. For the coefficients in the 
dipole approximation @ °f /i’ 9i we have gf: 

uiy = 6.34GeV, m A = 6.73GeV. 


The list of the processes and the corresponding decay 
widths is presented in the Table [T| 
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